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1. Introduzione 


La Facolta di Agraria dell’Universita di Palermo, sin dall’inizio del 2006, ha attivato una 
collaborazione con l’Associazione "Maison Shalom", in Burundi, nell’ambito di progetti inerenti lo 
sviluppo dell’agricoltura. L'Associazione dispone in Burundi, nelle province di Ruyigi e Cankuzo, 
di alcune aziende agricole che riescono a soddisfare soltanto il 30% del fabbisogno in alimenti della 
popolazione assistita dalla Maison Shalom. In questo quadro, la Facolta di Agraria ha avviato un 
programma di cooperazione, denominato "Progetto Burundi". L’impegno e stato finalizzato in 
primo luogo alia conoscenza del potenziale fabbisogno nutritivo e della frazione che attualmente 
viene soddisfatta. L’obiettivo finale della Facolta e quello di contribuire a realizzare forme di 
agricoltura innovative, che tengano conto anche della salvaguardia ambientale, nell’ottica di un 
progetto integrato finalizzato ad innescare processi di sviluppo autonomo, duraturo, sostenibile ed 
equo nel territorio. 

La Facolta di Agraria ha, inoltre, avviato uno studio volto ad ampliare le conoscenze sulle 
caratteristiche pedo-climatiche del paese, che ha portato nel 2007 alia stesura di un report di 
approfondimento sugli aspetti agro-meteorologi delle province di Ruyigi e Cankuzo (Bagarello et 
ah, 2007), e alia pubblicazione sulla rivista GEODERMA dell’articolo “Testing different 
approaches to characterize Burundian soils by the BEST procedure” (V. Bagarello, S. Di Prima, M. 
Iovino, G. Provenzano, A. Sgroi). A fine divulgativo e stato inoltre realizzato un sito internet 
(www.climateofburundi.org) , in cui sono consultabili e liberamente scaricabili i dati climatici di 
precipitazione e temperatura delle principali stazioni tenno-pluviometriche del Burundi. 

La presente tesi si inquadra nell’ambito delle attivita di ricerca svolte in Burundi, con il nuovo 
intento di analizzare un aspetto quanto mai emblematico dei nostri tempi. E’ stato affrontato uno 
studio sulle tendenze climatiche di lungo periodo riguardante il Burundi. L’analisi e stata svolta a 
scala nazionale sulla base della raccolta delle serie storiche di precipitazione e temperatura di piu 
lunga osservazione. 

Negli ultimi anni la percezione di attraversare una fase meteoclimatica anomala e aumentata 
trovando fondamento in studi realizzati a livello “planetario”. Gli studi globali rimandano tuttavia 
alia necessita di approfondimenti locali svolti sulla base di serie storiche disponibili. L’attivita di 
monitoraggio del clima non piu solo a livello globale, ma anche a scala locale, cerca si soddisfare la 
crescente domanda di informazioni relative alle ricadute che eventuali trend climatici avrebbero a 
livello locale sugli ecosistemi naturali ed agrari, e sull’ambiente in genere. 

Lo studio climatico del Burundi e stato sviluppato secondo i seguenti punti: 
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1 . Raccolta dei dati termo-pluviometrici; 

2. Selezione delle serie storiche (termo-pluviometriche) burundesi piu consistenti e 
possibilmente complete (cioe esenti da “ dati mancanti ” ); 

3. Ricostruzione dei dati mancanti; 

4. Individuazioni di eventuali trend a scala nazionale. 
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2. Raccolta delle serie storiche termo-pluviometriche burundesi 


Per condurre questa ricerca, sono stati acquisiti, presso l’lstituto Geografico del Burundi 
(IGEBU), i dati giornalieri di pioggia e temperatura dell’ aria registrati alle stazioni termo- 
pluviometriche di Bujumbura, Cankuzo, Gisozi, Gitega, Imbo, Karuzi, Kinyinya, Kirundo, 
Makamba, Mparambo, Mpota Tora, Muriza, Musasa, Muyaga, Muyinga, Nyanza Lac, Ruvyironza, 
Rwegura e Teza, per un totale di 20 stazioni termo-pluviometriche (Figura 1). In tabella 1 si riporta 
l’elenco complete delle stazioni termo-pluviometriche, e per ciascuna l’altitudine sul livello del 
mare, le coordinate geografiche, la precipitazione e la temperatura media dell’intero periodo di 
osservazione. 


Tabella 1 - Elenco delle principali stazioni termo-pluviometriche del Burundi. 


N. 

STAZIONE 

LAT. 

SUD 

LONG. 

EST 

ALTITUDINE 
[m s.l.m.] 

H a , mo [mm] 

T I°C1 

anno l v - x J 

1 

BUJ 

3°19' 

29° 19' 

783 

830.1 

24.0 

2 

CANK 

3°17' 

30°23' 

1652 

1215.1 

19.8 

3 

GIS 

3°34' 

29°4p 

2097 

1473.9 

16.6 

4 

GIT 

3°25' 

29°55' 

1645 

1208.1 

19.5 

5 

IMB 

3°10' 

29°20' 

834 

794.5 

24.1 

6 

KAR 

3°06' 

30°10' 

1600 

1158.8 

19.6 

7 

KIN 

3°39' 

30°20' 

1308 

1190.2 

21.9 

8 

KIR 

2°35' 

30°07' 

1449 

1071.1 

20.7 

9 

MAK 

4° T 

29°48' 

1424 

1269.9 

- 

10 

MPA 

2°49' 

29° 4' 

897 

925.6 

24.2 

11 

MPO 

3°44' 

29°34' 

2160 

1518.8 

15.4 

12 

MUR 

3°32' 

30°05' 

1616 

1156.5 

18.6 

13 

MUS 

4°00' 

30°06' 

1260 

1160 

21.6 

14 

MUYA 

3°13' 

30°32' 

1600 

1180.1 

- 

15 

MUYI 

2°51' 

30°07' 

1756 

1099.7 

19.8 

16 

NY AM 

2°52' 

30°E 

1651 

1361.9 

18.9 

17 

NY AN 

4°21' 

29°36' 

792 

1271.2 

23.4 

18 

RUV 

3°49' 

29°46' 

1822 

1317.4 

17.1 

19 

RWE 

2°54' 

29°30 

2226 

1678.2 

14.9 

20 

TEZ 

3°10' 

29°33' 

2440 

1618.8 

16.0 


II periodo di osservazione varia per ogni stazione e per entrambe le variabili meteorologiche; le 
osservazioni pluviometriche sono comprese tra il 1931 e il 2008, con un periodo massimo di 
osservazione pari a 74 anni solo per poche stazioni, le osservazioni tennometriche, invece, sono 
comprese tra gli inizi degli anni ’50 e il 2007, con un piu lungo periodo di osservazione pari a 35 
anni. Inoltre quasi tutte le serie presentano dati mancanti, riferiti all’intero anno o a singoli mesi. 
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Molte stazioni non presentano un numero minimo di anni di osservazione sufficiente a 
giustificare le analisi dei trend, condotte alia scala temporale mensile. Per questo motivo, le serie 
che non presentano almeno 15 anni di osservazione sono state escluse dal nostro database e non 
considerate nell’analisi. 



Figura 1 - Mappa della rete ternio-pluviometrica dell'Istituto Geogralico del Burundi (IGEBU). 

3. Selezione delle stazioni e archivio definitivo dei dati 

La scelta delle stazioni considerate, sulla base di alcune esigenze proprie dell ’analisi delle 
tendenze climatiche e stata effettuata sulla base dei seguenti criteri: 

• maggiore numero di anni di osservazione, 

• continuity dell’informazione termo-pluviometrica, 

• contemporaneity deH’informazione pluviometrica e tennometrica tra piu stazioni. 

In base a tali criteri sono state scelte le serie pluviometriche di 18 stazioni su 20, ovvero tutte le 
stazioni fatta eccezione per Musasa (MUS) e Nyanza Lac (NYAN). Con riferimento alle 
temperature, invece, sono state considerate 15 stazioni su 20. Le stazioni trascurate sono: Karuzi 
(KAR), Makamba (MAK), Muyaga (MUYA), Nyamuswaga (NY AM) e Nyanza Lac (NYAN). 

A partire dai dati giomalieri di ciascuna serie considerata, per ogni mese dell’intero periodo di 
osservazione, sono stati determinati l’altezza di pioggia mensile, espressa in millimetri, la 
temperatura media, massima e minima mensile dell’ aria, espressa in gradi centigradi e l’escursione 
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termica, pari alia differenza tra le temperature medie massime e minime. L’altezza di pioggia 
mensile e stata ottenuta dalla somma delle altezze di pioggia giornaliera. Le variabili di 
temperatura, invece, sono state calcolate come media delle temperature massime, minime e medie 
giomaliere. 


Le serie mensili di precipitazione e temperatura cosi ottenute sono state in seguito aggregate alle 
scale temporali stagionali. Per tale aggregazione temporale l’anno solare e stato suddiviso in due 
stagioni climatiche, una stagione umida compresa tra il primo ottobre e il 30 aprile (01/10 - 30/04) 
ed una stagione secca dallo 01/05 al 30/09. Di seguito, la figura 2 riporta, per le 18 stazioni 
considerate, l’andamento delle precipitazioni durante l’anno, da cui appare evidente la quasi assenza 
di precipitazioni durante la stagione secca (da maggio ad ottobre). 



Figura 2 - Andamento delle precipitazioni nel eorso dell’anno. 

In definitiva, quindi, sono state ottenute le serie di precipitazione, H (mm), relative al periodo 
umido, secco ed annuo, ottenute come somma dei volumi caduti nel periodo considerato. Alio 
stesso rnodo, le serie di temperatura massima, T max (°C), minima, T mm (°C), e media, T m .(° C), della 
stagione umida, secca ed annua, sono state ottenute come media dei valori mensili nel periodo di 
tempo considerato. L’escursione termica, AT (°C), annua e delle due stagioni e stata dedotta dalla 
differenza tra temperatura massima e minima annua e stagionale. 


4. Ricostruzione dei dati mancanti 


Seppure le serie termo-pluviometriche considerate sono tra le piu lunghe e continue, alcune di 
esse presentano dei periodi di mancata osservazione che interrompono la serie, rendendo vana ogni 
ulteriore analisi, per questo motivo si e ritenuto opportuno adottare una procedura per la 
ricostruzione dei dati mancanti. Tale ricostruzione e stata effettuata a scala stagionale, ritenendo 
poco attendibili i procedimenti di ricostruzione alia scala giornaliera. 
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II procedimento di ricostruzione dei dati mancanti e basato sul metodo dell’imputazione 
multipla, un approccio predittivo di analisi multivariata. 

Con il metodo dell’imputazione multipla (Rubin 1976, 1987), i valori mancanti di ciascuna 
variabile vengono ricostruiti utilizzando i valori delle altri variabili. In pratica, ogni singolo valore 
mancante viene sostituito da un set di possibili valori che rappresentano l’incertezza circa l’esatto 
valore da imputare. 

L’imputazione multipla consiste di una procedura sorprendentemente intuitiva, che consta delle 
seguenti fasi: 

1. i dati mancanti sono imputati in volte (generalmente da 3 a 5), generando cosi m dataset 

completi; 

2. gli m dataset completi sono analizzati usando i metodi dell’analisi statistica standard. 

3. i risultati dell’analisi sugli in dataset sono combinati per produrre il dataset completo. 

L’imputazione multipla aumenta l’efficienza della stima dei valori mancanti, in quanto la 
procedura di analisi dei diversi dataset imputati minimi zz a l’errore standard della ricostruzione. 
L’imputazione multipla permette non solo di mantenere la naturale variabilita dei dati mancanti, ma 
anche di includere l’incertezza causata dalla stima degli stessi. I valori imputati prodotti dal modello 
di imputazione, quindi, non van no intesi come valori "estratti", cosi come potrebbe essere inteso un 
singolo valore mancante, piuttosto, questo modello pennette di creare un data-set imputato che 
mantenga la variabilita complessiva della popolazione, e al contempo rispetti le relazioni con le 
altre variabili. 

Il metodo ha mostrato di produrre parametri di stima corretti che riflettono l’incertezza associata 
con la stima dei valori mancanti. Inoltre, l’imputazione multipla si e rivelato un metodo robusto e 
fomisce adeguati risultati anche in presenza di campioni di ridotte dimensioni o con alte percentuali 
di valori mancanti (Wayman, 2003). 

L’applicazione di tale metodologia alle serie di precipitazione e temperatura ha presupposto 
un’ analisi preliminare dei dati, volta a vcrilicare la percentuale di dati mancanti, il numero di anni 
di osservazione e nei casi possibili la disposizione geografica delle stazioni. In particolare per le 
serie di precipitazione si e scelto di non considerare come dataset originale l’insieme di tutte le 
stazioni, in quanto esse sono molto diverse in relazione soprattutto al numero di anni di 
osservazione. Considerarle insieme avrebbe determinate un’elevata percentuale di dati mancanti e 
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quindi avrebbe introdotto un’ulteriore incertezza nella ricostruzione, per niente funzionale ai Uni 
della nostra analisi. La ricostruzione, infatti, ha riguardato solo i dati mancanti compresi tra l’inizio 
e la fine del periodo di osservazione e non era volta all’estensione del periodo di osservazione. 

Per questi motivi sono stati individuati quattro gruppi di stazioni pluviometriche, le cui serie 
erano omogenee in tennini di periodo di osservazione e percentuale di dati mancanti, essi sono cosi 
costituiti: 

• Bujumbura, Gisozi, Makamba, Muyinga; 

• Karuzi, Mparambo, Muyaga, Bujumbura; 

• Gitega, Kinyinya, Mpota Tora, Muyaga, Ruvyironza, Rwegura, Teza, Karuzi; 

• Cankuzo, Imbo, Kirundo, Nyamuswaga, Karuzi, Muyaga. 

Alcune stazioni sono state inserite in piu gruppi al fine di non ottenere gruppi poco numerosi e 
per le loro caratteristiche di completezza, che li rendono dei buoni predittori nel processo di 
imputazione. 

Le stazioni termometriche sono state considerate, invece, in un unico gruppo, in relazione alia 
loro omogeneita nei confronti dei parametri sopra considerati, relativi alia percentuale di valori 
mancanti e al numero di anni di osservazione, abbastanza simile per tutte le stazioni considerate. 

Le serie climatiche sono state ricostruite mediante il pacchetto SPSS Missing Values 19, del 
software statistico SPSS, distribuito da IBM, scelta giustificata dalla facilita di implementazione del 
metodo e dalla agevole gestione dei dati di input ed output del software, anche se onnai, molti 
programmi statistici dispongono del pacchetto dell ’imputazione multipla. 

Le serie pluviometriche ricostruite presentano un periodo massimo di osservazione compreso tra 
il 1930 e il 2007 (periodo medio di osservazione pari a 52 anni - Figura 3, 4). Con 77 anni di 
osservazione la stazione di Bujumbura risulta essere la piu importante. Le serie termometriche 
presentano un periodo massimo di osservazione che va dal 1964 al 2007 e un periodo medio di 
osservazione pari a 33 anni (Figura 5, 6). 
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Figura 3 - Precipitazione nella stagione umida, numero di anni di osservazione prima e dopo la ricostruzione. 



Figura 4 - Precipitazione nella stagione secca, numero di anni di osservazione prima e dopo la ricostruzione. 
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Figura 5 - Numero di anni di osservazioni termometriche nella stagione umida prima e dopo la ricostruzione. 
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Figura 6 - Numero di anni di osservazioni termometriche nella stagione secca prima e dopo la ricostruzione. 
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5. Valutazione dei trend di precipitazione e temperatura a scala regionale 

Lo studio dei trend e stato condotto a scala regionale. In primo luogo e stato applicato il test non 
parametrico di Mann-Kendall a ciascuna serie delle grandezze climatiche prima elencate e 
successivamente e stato valutato il campo di signi flcati vita del trend a scala regionale mediante lo 
statistico medio regionale RAMK (S.Yue, C. Y. Wang, 2002). 

Lo statistico medio regionale RAMK (Regional Average Mann-Kendall), rappresenta le proprieta 
regionali del trend in quanto tiene conto dell’intero set di dati osservati nell’area di indagine. 
Inoltre, per tenere conto della struttura spaziale del trend in una specilica regione climatica e stato 
necessario valutame i campi di signi flcati vita, che risulta influenzato dalla correlazione seriale e 
dalla cross-correlazione tra i siti in una specifica regione. E’ noto che la presenza della 
autocorrelazione temporale {serial correlation ) e della cross-correlazione spaziale, in una serie 
climatica, influisce sulla signi licati vita dei test. Kulkami e von Storch (1995) hanno dimostrato che 
la presenza di autocorrelazione positiva comporta una piu alta probabilita di rigettare l’ipotesi nulla 
di “no trend” quando, invece, dovrebbe essere accettata. L’esistenza della cross-correlazione tra i 
siti conduce, invece, ad una ridondanza delle informazioni contenute in ogni sito; ad esempio, se 
due stazioni sono altamente correlate, la seconda stazione duplica meramente la prima e non 
fomisce maggiori informazioni. La cross-correlazione riduce l’effettivo numeri di siti. Cosi come 
l’auto correlazione, la presenza di cross-correlazione positiva, puo aumentare la probabilita di 
rigettare l’ipotesi nulla (assenza di trend) quando, invece, si dovrebbe accettare (Douglas, et al. 
2000). 

In particolare, nel test di MK, l’esistenza di correlazione seriale e di cross-correlazione altera la 
varianza della distribuzione nulla dello statistico di MK, ma non cambia la sua distribuzione 
normale (S.Yue, C. Y. Wang, 2002). Alla luce di cio, e stato sviluppato lo statistico regionale di 
Mann-Kendall, RAMK, la cui distribuzione nulla, si basa sul fatto che la distribuzione di probabilita 
congiunta di m variabili indipendenti nonnali e anch’essa distribuita nonnalmente, mentre, la 
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varianza del test statistico viene modificata o corretta mediante un fattore di correzione, che tiene 
conto deU’influenza della correlazione seriale e della cross-correlazione. 


Preliminarmente aU’applicazione del test statistico regionale, un’analisi delle serie storiche ha 
escluso la presenza di correlazione seriale significativa e quindi nessuna specifica e ulteriore 
correzione e stata prevista a riguardo. 

II test statistico regionale MK si sviluppa a partire dal test statistico non parametrico MK relativo 
ad ogni singola stazione, che e stato applicato ad ogni stazione ed ad ogni variabile considerata. 

II test non parametrico di MK si basa sullo statistico, S, dell ni to come segue: 

n—l n 

s = 2 2 sng(X 0 - XJ ( 1 ) 

i=l j=i + 1 

Dove Xk e il k-esimo valore della seria, n e la lunghezza della seria e la funzione sng(O) vale: 


( 1 per 0 > 0 

sng(0 ) = I 0 per 0 = 0 (2) 

(—1 per 0 < 0 

Lo statistico S, pertanto, rappresenta la differenza tra il numero di differenze positive e il numero 
di differenze negative nella serie data. Se e vera l’ipotesi nulla (Ho), cioe che non vi sia trend, ogni 
riordino del di dati originali ha la stessa probability; sotto tale ipotesi, Mann (1945) e Kendall 
(1975), dimostrarono che lo statistico S (per n > 8), e approssimativamente distribuito con legge 
nonnale, con media e varianza, dati rispettivamente da: 


E[S] = 0 

n(n — 1)(2 n + 5) 
Var(S ) = £ 


( 3 ) 

( 4 ) 


Per N>10, si puo poi dimostrare che la distribuzione dello statistico S’=S-1 approssima meglio la 
legge nonnale che non S; assumendo quindi che S’ sia distribuito nonnalmente, lo statistico 
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standardizzato di Mann-Kendall, Z, con media pari a zero e varianza unitaria, e computato come 
segue: 


Z = 


r S- 1 

—^=^= perS > 0 
y/var(S) 

< 0 per S — 0 


5 + 1 

perS < 0 

K yjVar (5) 


( 5 ) 


Un valore positivo di Z indica un trend crescente, mentre un valore negativo un trend 
decrescente. L’ipotesi H 0 (di assenza di trend) sara rifiutata se la probability p associata a Z e 
minore del livello di significativita del test, a. 


II p-value (minimo livello di significativita per il quale si rifiuterebbe Ho) dallo statistico 5 di 
MK del campione di dati, puo essere stimato usando la funzione di distribuzione cumulata normale: 


P = 0.5 — d>(|Z|) (6) 

Z 

lft 2 

Q(\Z\) = —J e 2 dt (7) 


o 

S q p < a , il segno di Z definisce l’andamento del trend, crescente per 5 positivo, decrescente in 
caso contrario. 


Per valutare il campo di significativita del trend a scala regionale, e stato adottato lo statistico 
medio regionale MK ( RAMK ) proposto da Yue & Wang (2002), la cui distribuzione nulla e derivata 
sulla base della distribuzione congiunta di m variabili indipendenti. 


La somma dello statistico MK relativo a m siti indipendenti e data da: 


m 


$SUM = 




Dove Sk e lo statistico MK relativo a ciascun sito k. 


( 8 ) 
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La distribuzione nulla dello statistico Ssum e normale; infatti, la distribuzione congiunta della 
somma di m variabili indipendenti distribuite con legge normale, e anch’essa normale {Hogg & 
Craig, 1978 in Yue & Wang, 2002): 


m 

J ($sum) = Urw 

k = 1 

( m m \ 

/ Mfc - / a k I 

k=l k=l / 


( 9 ) 


( 10 ) 


Dove jUk e sono rispettivamente la media e la varianza dello statistico MK del siti /ceil 

2 

simbolo ~ indica eguaglianza asintotica alia distribuzione normale, N. In pratica //a e <r* sono posti 
equivalenti rispettivamente alia media E[Sk ] e alia varianza Var(Sk) campionaria. 


Se Ssum e diviso per la costante m si ottiene lo statistico RAMK S=Ssum hn anch’esso distribuito 
normalmente: 



La media e la varianza campionaria dello statistico RAMK S valgono rispettivamente: 
E[S\ = 0 

m 

var (s) =-^Y u var( - Sk ) 

k = 1 

Lo statistico RAMK standardizzato S risulta quindi il seguente: 


(ID 


(12) 

(13) 


5-1 


per S > 0 


Z = < 


Var(S ) 

0 per S — 0 


5 + 1 


per 5 < 0 


Var(S ) 


(14) 
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La distribuzione nulla dello statistico RAMK deriva da Yue & Wang (2002), e la stessa di quella 
dedotta da Douglas et al. (2000) in forza del teorema limite centrale. Un pregio della presente 
derivazione e quello che, essendo la normalita della distribuzione nulla derivata da una 
distribuzione congiunta, vengono meno i limiti imposti al numero delle stazioni dal teorema limite 
centrale, per cui la normalita e, a rigore, asintotica. La distribuzione nulla dello statistico RAMK puo 
essere applicata ad una rete che a persino due stazioni di osservazione (m=2). 

Come gia precedentemente accennato, l’esistenza di cross-correlazione tra le stazioni riduce 
l’effettiva dimensione del campione e accresce la varianza dello statistico RAMK (Douglas et al. 
2000), lo statistico regionale Z e la varianza di S ,Var(S), sono stati corretti. 

La varianza di S in presenza di cross-correlazione nella serie e modificata nel modo seguente: 
Lar(s) = Lar(5)r7 c (r] c = 1 + (m — l)p fcfe+1 ) (15) 


In cui if e il fattore di correzione per la cross-correlazione ep k k+1 e la media dei coefficienti di 
cross-correlazione di tutti i siti considerati. Esso presenta la seguente espressione: 


Pk,k+ 1 


2 £ 


m—l 

K=l 


srm—k _c 

2.1=1 Pk,k + 1 


m(m — 1) 


(16) 


Dove p c k k+1 e il coefficiente di cross correlazione tra il sito k e il sito k+1. 
Lo statistico standardizzato RAMK, Z r , e, infine pari a: 


~ z '-j+ (17) 

Infine, supposta disponibile una serie di dati osservati, X (=x / , X 2 , ...,x„), l’intensita del trend, 
bo, puo essere computata usando l’approccio di Theil (1950) e Sen (1968), TSA, la cui espressione 
e la seguente: 

b 0 — Median p- j V/ < j (18) 

Dove xi e 1’1-esimo valore del campione di datiX. 
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6. Risultati 


Le Tabelle 2 e 3 riportate di seguito espongono la sintesi dei risultati dell’analisi climatica della 
regione centro africana. In tabella 2 sono mostrati i valori regionali dello statistico medio regionale 
RAMK, per tutte le variabili indagate e le scale temporali considerate (stagione umida, StJUm, 
stagione secca, St_Sec, ed anno), in riferimento ad un livello di significativita del trend pari al 5%. 


Tabella 2 - Valori dello statistico medio regionale RAMK, Zr, ed il p-value (in grigio sono 
riportati i valori non significativi al 5%) 



St _Um 

StSec 

ANNO 

Zr 

p-value 

Zr 

p-value 

Zr 

p-value 

T 

-* max 

2.618 

0.009 

3.796 

0.0001 

3.309 

0.001 

T 

x min 

3.630 

0.0003 

3.837 

0.0001 

3.577 

0.0003 

T m 

2.679 

0.007 

4.146 

3E-05 

3.211 

0.001 

AT 

0.205 

0.837 

0.052 

0.958 

0.430 

0.667 

H 

-0.972 

0.331 

-1.245 

0.213 

-1.395 

0.163 


Le temperature massime, minime e medie mostrano una significativa tendenza complessiva 
aU’incremento, per tutte le scale temporali considerate. Al contrario, con riferimento alle serie di AT 
e di precipitazione, i risultati dell’analisi non hanno evidenziato, sia alia scala stagionale che a 
quella annua, un trend significativo per l’intero Burundi. 

In tabella 3, per ciascuna variabile climatica sono riassunti i valori regionali del trend computati 
con l’approccio di Theil applicato agli scarti rispetto alia media dei dati di tutte le stazioni 
considerate. 


Tabella 3 - Valori regionali dei trend: T, max , T mi „, T m e AT [°C anno" 1 ], H [mm anno' 1 ], (in 
grigio sono riportati i valori non significativi al 5%) 



St Um 

St Sec 

ANNO 

T 

-* max 

0.016 

0.023 

0.017 

T min 

0.015 

0.021 

0.017 

T m 

0.013 

0.021 

0.015 

AT 

-0.0003 

0.001 

0.001 

H 

-0.546 

-0.347 

-0.983 


Alla scala stagionale l’incremento della temperatura appare imputabile principalmente alia 
stagione secca, con una tendenza d’incremento di 0.02 °C/anno delle temperature medie della 
stagione. 
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L’assenza di trend per la variabile AT, puo essere giustificato dal fatto che le temperature 
minime e massime mostrano un trend positivo alle diverse scale considerate dello stesso ordine di 
grandezza come riportato in tabella 3. 

La figura 7, che riporta l’andamento temporale degli scarti dalla media delle precipitazioni, 
evidenzia l’assenza di un trend significative. 

Nelle figure 8, 9 e 10 e, invece, riportato l’andamento cronologico degli scarti delle temperature 
medie, massime e minime stagionali ed annue. I grafici mostrano la chiara tendenza aU’incremento 
delle temperature. 

7. Conclusioni 

In definitiva dall’analisi delle serie storiche di precipitazione, registrate in Burundi dal 1939 al 
2007 in 18 stazioni termo-pluviometriche, non emerge una variazione statisticamente rilevante, a 
scala regionale, dei valori di precipitazione, per tutte le scale considerate. Solo poche stazioni 
mostrano una tendenza alia diminuzione. Cio non implica che il global warming non causi in 
Burundi delle variazioni sui regimi di precipitazione, ma indica che i cambiamenti in atto non sono 
statisticamente rilevanti e/o rilevabili con i dati disponibili. 

L’analisi sulle serie storiche di temperatura, registrate dal 1964 al 2007, in 15 stazioni 
termometriche del Burundi mostra, invece, una tendenza aU’incremento delle temperature che, a 
scala annua, assume un valore pari a 0.015 C°/anno. 

Con riferimento alia precipitazione, le conclusioni dello studio sono concordi con altre indagini 
eseguite da altri autori nell’ Africa centrale: New et al., 2006, per l’ovest e il sud dell’Africa; 
Conway et al., 2008, per l’Africa sub-Sahariana; Aguilar et al., 2009, per l’Africa centrale, 
orientale, la Repubblica di Guinea e lo Zimbabwe; Mazvimavi, 2010, per lo Zimbabwe. 

Le dinamiche climatiche individuate, che comporterebbero una modifica del bilancio idrologico 
nella direzione dell’aumento delle perdite, necessita di un approfondimento volto a valutame gli 
effetti sugli ecosistemi naturali ed in particolare agrari, del Burundi. 
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Figura 7 - Andamenti degli scarti dalla media delle serie di precipitazione alle scale temporali umida, secca ed annua. 
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Figura 8 - Andamenti degli scarti dalla media delle serie di temperatura massima, minima e media nella stagione umida. 
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Figura 9 - Andamenti degli scarti dalla media delle serie di temperatura massima, minima e media nella stagione secca. 
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Figura 10 - Andamenti degli scarti dalla media delle serie di temperatura massima, minima e media annua. 
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La Facolta di Agraria dell'Universita di Palermo, sin dall'inizio del 2006, ha attivato una 
collaborazione con l'Associazione "Maison Shalom", in Burundi, nell'ambito di progetti inerenti lo 
sviluppo dell’agricoltura. L'Associazione dispone in Burundi, nelle province di Ruyigi e Cankuzo, di 
alcune aziende agricole che riescono a soddisfare soltanto il 30% del fabbisogno in alimenti della 
popolazione assistita dalla Maison Shalom. In questo quadro, la Facolta di Agraria ha avviato un 
programma di cooperazione, denominato "Progetto Burundi". L’impegno e stato finalizzato in primo 
luogo alia conoscenza del potenziale fabbisogno nutritivo e della frazione che attualmente viene 
soddisfatta. L’obiettivo finale della Facolta e quello di contribuire a realizzare forme di agricoltura 
innovative, che tengano conto anche della salvaguardia ambientale, nell’ottica di un progetto integrato 
finalizzato ad innescare processi di sviluppo autonomo, duraturo, sostenibile ed equo nel territorio. La 
Facolta di Agraria ha, inoltre, avviato uno studio volto ad ampliare le conoscenze sulle caratteristiche 
pedo-climatiche del paese, che ha portato nel 2007 alia stesura di un report di approfondimento sugli 
aspetti agro-meteorologi delle province di Ruyigi e Cankuzo (Bagarello et al., 2007), e alia 
pubblicazione sulla rivista GEODERMA dell’articolo “Testing different approaches to characterize 
Burundian soils by the BEST procedure” (V. Bagarello, S. Di Prima, M. Iovino, G. Provenzano, A. 
Sgroi). A fine divulgativo e stato inoltre realizzato un sito internet (www.climateofburundi.org), in cui 
sono consultabili e liberamente scaricabili i dati cbmatici di precipitazione e temperatura delle 
principali stazioni termo-pluviometriche del Burundi. La presente tesi si inquadra nell’ambito delle 
attivita di ricerca svolte in Burundi, con il nuovo intento di analizzare un aspetto quanto mai 
emblematico dei nostri tempi. E’ stato affrontato uno studio sulle tendenze climatiche di lungo periodo 
riguardante il Burundi. L’analisi e stata svolta a scala nazionale sulla base della raccolta delle serie 
storiche di precipitazione e temperatura di piu lunga osservazione. Negli ultimi anni la percezione di 
attraversare una fase meteoclimatica anomala e aumentata trovando fondamento in studi realizzati a 
livello “planetario”. Gli studi globab rimandano tuttavia alia necessity di approfondimenti locali svolti 
sulla base di serie storiche disponibili. L’attivita di monitoraggio del clima non piu solo a livello 
globale, ma anche a scala locale, cerca si soddisfare la crescente domanda di informazioni relative alle 
ricadute che eventuali trend climatici avrebbero a livello locale sugli ecosistemi naturali ed agrari, e 


sull’ambiente in genere. 



